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In  this  paper  two  new  series  of  search  designs  with  very  small 
number  of  treatments  are  presented  for  3™  factorial  experiments.  The 
first  series  of  designs  can  search  one  nonzero  two  factor  interaction 
and  estimate  it  along  with  the  general  mean  and  the  main  effects. 

The  second  series  of  designs  can  search  one  nonzero  three  factor 
interaction  and  estimate  it  along  with  the  two  factor  and  lower  order 
interactions. 


Short  Running  Title:  Search  Designs  for  3™  Factorials 

AMS  1970  Subject  Classifications:  Primary  and  Secondary:  62K15 

Keywords  and  Phrases:  Estimability,  Interactions,  Main  Effects, 
Resolution  III  and  V  Plans,  Search  Designs. 


*The  work  of  the  first  author  is  sponsored  by  the  Air  Force  Office  of  □ 

□ 


Scientific  Research  under  Grant  AF0SR-86-0048. 

**The  work  of  the  second  author  is  supported  by  a  fellowship  from  the 
University  of  California,  Riverside. 

|  — 

Dbt 


Codes 


I 


1,;  and/or 
Special 


1 .  Introduction 


Consider  the  linear  model 

(1)  E(£)  -  XA  jlj  +  X2  J2, 

(2)  V  (2.)  -  0*1, 

where  (Nxl)  is  a  vector  of  observations  and  for  i  ■  1,2,  (Nxv^) 

are  known  matrices,  B.(v  xl)  are  vectors  of  fixed  parameters  and  a2  is 

— i  i 

a  constant  which  may  or  may  not  be  known.  Moreover,  J1A  is  completely 

unknown,  but  we  have  partial  information  about  J5 2-  We  know  that  at 

most  p  elements  of  jhj  are  nonzero  and  the  remaining  elements  are 

negligible,  where  p  is  a  nonnegative  integer  which  may  or  may  not  be 

known.  In  this  paper  we  assume  p  is  known  to  be  1.  However,  we  do 

not  know  exactly  which  element  of  *8  nonzero.  The  problem  is  to 

search  the  nonzero  element  of  and  draw  inference  on  it  in  addition 

to  the  elements  of  .  Such  models  are  called  search  linear  models 

and  were  introduced  in  Srivastavs  (1975).  We  want  Xj  and  X2  to  be 

such  that  the  above  problem  can  be  resolved;  the  underlying  design 

corresponding  to  XA  and  X2  is  called  a  search  design. 

In  a  3™  factorial  experiment  the  treatments  are  denoted  by 

(a  ,...a  ),  a  *  0,1,2;  the  factorial  effects  are  deonoted  by 
1  ml 

C1  C 

F,  . .  .F  m  c.  ■  0, 1 ,2;  the  observation  corresponding  to  the  treatment 
1  ml 

(a  ,...,a  )  is  denoted  by  y  (a  ,...., a  ).  The  expectation  form  of  the 


model  is 


-3- 


where  the  values  of  depend  on  and  and  are  given  in  Table  1. 

Table  1 

The  values  of  b^ 


\ci 

ai \. 

0 

1 

2 

0 

1 

-1 

1 

1 

1 

0 

-2 

2 

1 

1 

1 

We  now  consider  the  following  two  situations: 

SI:  The  vector  consists  of  the  general  mean  and  the  main  effects 
and  the  vector  consists  of  the  2-factor  interactions.  The 
3-factor  and  higher  order  interactions  are  assumed  to  be  zero. 

We  assume  m  _>  3. 

S2:  The  vector  consists  of  the  general  mean,  the  main  effects  and 
the  2-factor  interactions  and  the  vector  consists  of  the 
3-factor  interactions.  The  4-factor  and  higher  order  inter¬ 
actions  are  assumed  to  be  zero.  We  assume  m  4. 

For  given  N  treatments,  we  can  write  the  equation  (3)  in  the  form  of 
the  equation  (1)  for  both  SI  and  S2.  Let  DA  be  a  design  with  (l+2m) 
treatments  as  the  treatment  with  all  factors  at  level  2,  treatments 
with  the  ith  (i  *  1 , . . . ,m)  factor  at  levels  of  0  and  1  and  the  other 
factors  at  level  2.  We  know  that  under  DA  we  can  estimate  the 


elements  of  8^  in  SI  with  the  assumpution  that  ■  _0  (i.e.,  a 
Resolution  III  plan).  Let  D2  be  a  design  with  treatments 

as  (l+2m)  treatments  in  ,  treatments  with  the  (i,j)th  factors 

—  ,m ,i  <  j)  at  levels  (0,0),  (0,1),  (1,0),  (1,1)  and  the  other 

factors  at  level  2.  We  know  that  under  D2  we  can  estimate  the 

elements  of  ^  in  S2  with  the  assumption  that  •  0^  (i.e.,  a 

Resolution  V  plan) .  Let  D3  be  a  design  with  m  treatments  as 

treatments  with  the  ith  (i“l,...,m)  factor  at  level  0  and  the  other 

factors  at  level  1.  We  prove  that  the  design  consisting  of 

(2) 

treatments  in  D|  and  D3  and  the  design  D  consisting  of  the 
treatments  in  D2  and  D3  are  in  fact  search  designs  in  SI  and  S2 , 
respectively. 

In  all  Taguchi  design  methods.  See  Taguchi  and  Wu  (1985),  popular 
in  statistical  quality  control  experimentations,  the  higher  order 
interactions  (2-factor  and  higher  order  in  most  plans)  are  assumed  to 
be  zero.  A  few  of  those  higher  order  interactions  may  have 
significant  effect  on  the  optimal  experimental  condition.  The  use  of 
search  designs  may  be  a  potential  tool  in  improving  upon  the  Taguchi 
design'  methods. 

There  is  a  vast  literature  available  in  the  construction  of 
search  designs  for  2™  factorial  experiments.  Near  minimal  resolution 
IV  plan  which  permit  search  and  estimation  of  three  or  fewer  nonzero 
two  factor  interactions  for  pm  factorial  experiments  are  available  in 
Anderson  and  Thomas  ( 1980)  under  the  assumption  that  3-factor  and 
higher  order  interactions  are  all  zero.  Chatter jee  and  Mukherjee 
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(1986)  presented  search  design  for  sr  x  */m  where  s  and  w  are  any 

positive  integers,  under  the  assumption  that  3-factor  and  higher  order 

(2) 

interactions  are  all  zero.  Our  Dv  in  S2  is  therefore  totally  new 
and  there  is  no  other  competitor  in  literature.  Our  design  in  SI 

is  although  new  but  has  competitors  in  Chatterjee  and  Mukherjee  (1986) 
and  in  Anderson  and  Thomas  (1980).  However,  the  design  has  an 

edge  over  designs  in  those  papers  in  terms  of  the  smaller  number  of 
treatments.  This  can  be  seen  by  taking  examples  4.5  and  4.6  in 
Chatterjee  and  Mukherjee  with  r  -  m  and  s  ■  3  and  comparing  with  the 
number  of  treatments  (l+3m)  in  .  In  SI ,  the  design  in  the 

example  4.5  in  Chatterjee  and  Mukherjee  has  ( l+6m)  treatments  with  3m 
more  treatments  than  in  ^ .  In  example  4.6  (m-3),  Chatterjee  and 
Mukherjee  has  8  more  treatments  than  in  Indeed,  both  Chatterjee 

and  Mukherjee,  Anderson  and  Thomas  designs  can  be  used  in  factorial 
experiments  other  than  3™. 

2.  Preliminary  Results 

We  first  introduce  the  following  notations. 

gy 

y^  «  The  observation  corresponding  to  the  treatment  with  levels  of 
all  factors  except  i  and  j  are  a,  the  level  of  the  ith  factor 
is  8  and  the  level  of  the  jth  factor  is  y. 

g 

y  *  The  observation  corresponding  to  the  treatment  with  levels  of 
all  factors  except  i  are  a  and  the  level  of  the  ith  factor  is 


8. 


-6- 


ya  “  The  observation  corresponding  to  the  treatment  with  levels  of 
all  factors  are  a. 

c  c 

S(c  «u)  ■  Sum  of  factorial  effects  F,  ...F  m  with  c»u. 
i  1  m  i 

c  c 

S(c *u,c *v)  -  Sum  of  factorial  effects  F.  ...F  °  with  c.“u  and  c  *v. 
i  j  1  m  i  J 

We  now  present  the  minimum  variance  unbiased  estimators  (MVUE)  of 


S(c  «u)  and  S(c  *u,c  -v)  under  (1)  and  (2)  with  B  -0  for  both  designs 

•  J  ^ 

Dj  and  D2 . 


Table  2 


MVUE' 8  of  S(ci«u)  and  S(c^-u,c  •“v)  for  D  and 

—  . .  — —  .  — .  — ,  .  V.  -  ,  -  ■ ,  ■  — 


sign 

Parameter 

MVUE 

2  S(c±-1) 

6  S(Ci-2) 

y2~y21 

y2"2y2i+y2i 

»2 

4  S(c1-l,cj»l) 

12  S(c1=l,cj»2) 

12  SCc^.Cj-l) 

36  S(c1-2,cj-2) 

y2  y2i  y2j  y2ij 

(y  _v0  )+(y0  _y00  )_2(y1  -?y°*  ) 

'y2  y2imy2j  y2ij '  ^y2j  y2ij' 

(y2-y5j)+(y5i-y5!j>-2(y21-yi;J) 

<y2ty2l'2y21>+(y2j"2y21j+y21j) 

~2 ( y2 j ”2y2i j  +y2i j  ^ 

The  requirement  on  Xj  and  X2  for  a  design  to  be  a  search  design 
(see,  Srivastava  (1975))  is  that  for  (^jmodels 

<4>  e(£)  -  +  x^  e*1*,  i-i . (v|), 

where  X2  (Nx2)  is  a  submatrix  of  X2  and  $2^  is  a  (2x1)  subvector 
of  J52,  the  parameters  and  are  unbiasedly  estimable.  Both 

°I  ^  and  consist  of  two  component  designs  namely  (DlfD3)  and 

(^2,1)3).  For  u*l ,2,3,  we  denote  the  observations  corresponding  to  D 
by  and  write  (4)  as 


e<v  -  >'uii1+41)-?<2)-1-1 . (v22)- 


For  u*l,2,  we  have 


Rank  X  ,  *  v 
ul  1 


E[-X31X^2u+i3]  -  [X32)-X31XuJXu2)l®21)- 


We  denote 


/(i)  -  Fx(i)  y  v"1Y<i)l  19-.,  fVz) 

*u  ^X32  X31XulXu2  1,2,1“1 . (  2-'* 


For  u**l,2,  the  requirement  for  the  design  to  be  a  search  design 


is  that 


Rank  W(i)  -  Rank  W(i)’w(i)  -  2,  i-1 , . . . ,(V?) . 
u  uu  ’  1  \  2 ' 


(i)  /fi 

It  16  to  be  noted  that  columns  of  W  ( rows  and  columns  of  W  w  ‘ 


correspond  to  a  pair  of  elements  of  Since  both  designs  D  and 


are  balanced  arrays  of  full  strength,  it  can  be  seen  that  for 


two  pairs  of  elements  of  g,  isomorphic  under  the  symmetric  group  of 


permutations  of  degree  m  on  the  set  {l,2,...,m},  the  W^  matrices  are 


identical  except  for  permutations  of  rows  and  columns.  We  now  present 
in  Table  3  and  4  the  nonisomorphic  pairs  of  elements  in  for  SI  and 


S2.  We  note  that  2-factor  interactions  can  be  any  of  3  types  F^F  * 


F^Fj ,  F2F2  for  different  values  of  i,j  without  assuming  the  restric¬ 
tion  i  <  j.  Similarly  3-factor  interactions  can  be  any  of  4  types 

F^F^F,  ,  F2F,F.  ,  F2F2F,  ,  F2F2F2  for  different  values  of  i,j,k  without 
i  j  k  i  j  k  i  j  k  i  j  k 

assuming  the  restriction  i  <  j  <  k. 


Table  3 

The  nonisomorphic  pairs  of  elements  in  for  SI 


(FiFJ'FlFo> 

<FiFrFiFu) 

(F2Fj,F2F 
i  j  u  v 

(FiFrFuFv> 

<FiFrFuF^ 

<FiFj  ,FiFj 

<FiFrFiFj> 

(FiFj •FiF„) 

(FSFj -Fft 

<FiFj •F1FU) 

<FiFj  ,FjFi> 

(FiFj  ,FjFu 

<FiFJ •FUF1> 

<Fi2FJ>FJFu> 

(fJfj.fV: 

(FiFj'FuFv) 

<FiFj ■FuFi) 

(FiFrFiFu 

<FlYriFj> 

(F2F  ,F2F  ) 

i  J  “] 

(f2f2  f2f2: 

i  j  U  V 

ess 


ran*r*'^TT^"v^"»'wu  w  u 
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Table  4 

The  nonisomorphic  pairs  of  elements  in  for  S2 


It  follows  from  (7)  that  w(i)g(i)  fl  8et  Qf  two  equations 

u  u  — 2 

In  two  elements  of  and  the  checking  of  Rank  ■  2  can  be 

—2  u  u 

done  by  showing  two  independent  unbiasedly  estimable  equations  in 
elements  of  0^*^ .  We  approach  this  problem  by  considering  Table  2, 
the  equation  (5)  for  u  -  3  and  the  equation  (5)  for  the  treatment  with 
all  factors  at  level  2. 


3 .  Main  Results 

We  now  present  our  two  main  results  and  their  proofs. 

Theorem  1 .  The  design  ^  is  a  search  design  for  SI. 

Proof .  The  proof  consists  of  showing  two  independent  unbiasedly  esti¬ 
mable  parametric  functions  of  the  elements  in  every  pair  in  Table  3. 

We  explain  the  nature  of  the  proof  by  considering  only  one  out  of  21 

pairs  in  Table  3  for  the  lack  of  space.  We  consider  the  model  (4) 

(1)  2  2 

with  the  elements  of  as  (F^Fj,F^Fj).  From  the  design  and  the 

parametric  function  S(c  **1)  in  Table  2,  it  follows  that  the  parametric 


u 


functions  (i)F  +F  F  ,  (ii)  F  +F  F,  and  (iii)  F  ,  u*i,j,  are  unbiasedly 

1  1  j  j  1  j  u 

estimable.  Again,  from  the  design  D,  and  S(c  *2)  in  Table  2,  it  can  be 

1  u 

seen  that  (iv)  F^+F^ Fj  ,( v)F2+F2 F2  and  (vi)  F  2 ,  u*i,  j,  are  unbiasedly 
estimable.  For  the  treatments  i,j  and  the  other  treatments  in  ,  we 


find  from  the  equation  (5)  that  (vii)  U-F1+F^-2Fj-2F^FJ2 ,  (viii)  y-Fj 
+F2-2F2-2j*’2 Fj  and  (ix)  y-2(F2+F2)  are  unbiasedly  estimable.  For  the 


LdJM. 


A.n.lUV 


•.vVQ 


-11- 


treatment  with  all  factors  at  level  2,  we  find  from  (5)  that  (x)  p+F^ 
+Fj+F^+Fj+Fi Fj+F^Fj  unb^-asedly  estimable.  From  (ix)-(vii)+3(iv)- 
(1),  it  follows  that  -F^Fj+5F^F^  is  unbiasedly  estimable.  From  (vii) 
+(viii)+(x)-3(ix)-6((iv)+(v)) ,  we  find  that  F^  -15F*F*  is  unbiasedly 
estimable.  We  thus  have  displayed  two  independent  unbiasedly  esti¬ 
mable  parametric  functions  in  F^F^  and  F2F  ^  .  The  checkings  for 
other  pairs  in  Table  3  are  similar  to  the  above.  This  completes  the 
proof  of  the  theorem. 

(2) 

Theorem  2.  The  design  D  is  a  search  design  for  S2 . 

Proof .  The  proof  consists  of  showing  two  independent  unbiasedly 


estimable  parametric  functions  of  the  elements  in  each  of  68  pairs  in 


Table  4  and  we  explain  the  nature  of  the  proof  by  considering  the  only 

pair  (F^F^  F^  .F^F^  F^) .  From  the  design  and  the  parametric  functions 

S(cv=l,  cw°l)  it  follows  that  (i)  F^F^  with  (v,w)  *  ( j ,k)  and  (ii) 

2 

FiFk+FiFjFk  are  unbiasedly  estimable.  From  the  parametric  functions 
S(cv  =  l,  c^'2) ,  S(cv*=2,  cw=l)  and  S(c^**2,  0^=2),  it  follows  that  (iii) 
F*F  with  (v-w)  *  Ci.J).  U,u),  (j,u),  (iv)  F2F+F2F  F,  ,  (v) 

V  W  i  j  i  j  k 

pjF-.+pfV..  (vi)  F2F  +F2F2F  ,  (vii)  F2F  +F?F?F  ,  (ix)  F2F2  with 
ikijk  iuiju  juiju  vw 

(v,w)  *  (i,j)  and  (x)  F2F2+F2F2F  are  all  unbiasedly  estimable.  It 

*  J  *  J  * 

follows  from  S(c  =1)  and  S(c  =2)  that  (xi)  F  ,  v  *  j.k.u,  (xii) 
v  v  v 

VF^+FA+FiFA'  <xiil)  V^WrtVk-  (xiv)  Vfi2V 


FJ  Fu+Fi  Fj  Fu  ’  <KV)  Fv’  v  *  i’^’  <*Fl>  Fi‘+Fi‘rj+Fl‘Fk+Fl‘F„+Fi‘Fj+Fl'FjFk+ 

F2F2F  ,  (xvii)  F2+F2F  +F2F2+F2F2F  are  unbiasedly  estimable.  For 
iju  jjuijiju 


2  2 
1  T?  T 


,2„2.r,21 


V. 


111  Mu  »~i  i 


the  treatments  i,j,k,u  and  the  other  (for  m  >  4)  treatments  in  D3  ,  we 

get  from  (5).  (ix),  (xi)  and  (xv)  that  (xviii)  p+F2-2F2-2F2F2 ,  (xix) 

'  i  j  i  J 

M_Fj+Fj"2Fi+2FiFj-2?2FJ’  <XX)  li-Fk-2F2-2F2+2F2Fk+4F2F2,  (xxi)  p-Fu 

-2F2-2F2+2F?F  +2F2F  +4F2F2-4F2F2F  and  (xxii)  p-2F2-2F2+4F2F2  are  all 
i  j  iu  ju  ij  iju  M  i  j  ij 

unblasedly  estimable.  (We  infact  will  not  use  the  equation  (xxii).) 

From  (xix)  -  (xviii)  +3  (vii)  -3  (vi)  -3  (v)  -6  (iv)  -  (ii)  -3  (xvii) 

+3  (xvi)  +  (xii)  and  (xx)  -  (xviii)  -9  (x)  -3  (vi)  -6  (v)  -3  (iv)  - 

(ii)  +3  (xvi)  +  (xiii)  we  find  that  -6F^ F.  F,  and  -6F2F,  F.  -9F2F2F 

ijk  i  j  k  iju 

are  unblasedly  estimable.  This  displays  two  independent  linear 

functions  of  F2F2F.  and  F2F2F  which  are  unblasedly  estimable.  The 
ijk  Iju 

checkings  for  the  other  pairs  in  Table  4  can  be  done  similarly.  This 
completes  the  proof  of  the  theorem. 


4.  Concluding  Remarks. 

a.  The  design  has  ( l+3m)  treatments  and  the  design 

2 

has  ( l-Hn+2m  )  treatments.  Minimal  resolutions  III,  V  and  VII 
plans  require  (l+2m),  (l+2m2)  and  (l+2m2+8(™))  treatments. 

b.  The  design  ^  can  search  one  nonzero  two  factor  inter¬ 
action.  A  natural  question  comes  up,  "can  search  one 

nonzero  3-factor  or  higher  order  interaction?”  The  answer  is 

"NO".  For  example,  in  case  m*»5  if  we  consider  the  model  (4) 

(i)  2  2 

with  the  elements  of  J8£  as  (FAF2  ,F3Fi,F5) ,  we  can  not  find 


two  independent  unblasedly  estimable  linear  functions  of  F1F2 
and  F3F$F5.  There  are  in  fact  many  such  pairs. 
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By  calling  the  level  2  as  the  level  0  (the  method  of 
collapsing  levels),  and  omitting  the  replicated  treatments, 
we  get  essentially  the  two  series  of  search  deisgs  obtained 
in  Srlvastava  and  Ghosh  (1976),  Srivastava  and  Gupta  (1979) 
for  2m  factorial  experiments.  However,  the  designs  thus 
obtained  have  more  strength  in  the  sense  that  they  can  search 
any  nonzero  i-f actor  or  higher  order  interactions,  where 
i  *  2  and  3. 

This  research  started  from  an  unpublished  technical  report  of 
Ghosh  (1985)  and  the  motivation  was  to  reduce  the  number  of 
treatments.  The  designs  ^  and  show  the  success  in 

our  effort. 
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